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Three-nucleon scattering equations with irreducible three-nucleon force are solved in momentum- 
space. The Coulomb interaction between the two protons is included using the method of screening 
and renormalization. The need for the renormalization of the scattering amplitudes is demonstrated 
numerically. The Coulomb and Urbana IX three-nucleon force effects on the observables of elastic 
proton-deuteron scattering and breakup are studied. 
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I. INTRODUCTION 

Nucleon-deuteron scattering has been extensively stud- 
ied in the past with the aim of testing various nuclear in- 
teraction models, i.e., nucleon-nucleon (NN) potentials 
and also three-nucleon (3N) forces. While the theoretical 
description of neutron- deuteron (nd) scattering is con- 
siderably simpler and a number of highly sophisticated 
calculations exists with realistic NN potentials and 
3iV forces, the experimental data are much more abun- 
dant and precise for proton-deuteron (pd) reactions. In 
that case there are only few realistic calculations above 
the deuteron breakup threshold. Configuration-space 
treatment based on the Kohn variational principle [a, 7| 
uses the local Argonne Vis (AV18) NN potential |8| 
and the irreducible Urbana IX (UIX) 3N force @ but 
is limited so far to elastic scattering at energies below 
50 MeV in the center-of-mass (cm.) system; only very 
recent work [To| solving integral Faddeev equations in 
configuration-space provides pd breakup results at pro- 
ton lab energy 13 MeV. In contrast, the momentum-space 
treatment pi], KL2| based on the screening and renormal- 
ization method PJ] has lead to pd elastic scattering 
and breakup results for energies up to pion-production 
threshold; this method can be used with nonlocal NN 
potentials as well, but so far only an effective 3N force 
due to the virtual excitation of a nucleon to a A isobar 
3 has been included. The aim of the present work is to 
overcome that limitation, i.e., to extend the technique of 
Refs. [ill, E3 to include also an irreducible 3N force. 

In Sec. [II] we derive three-particle scattering equations 
including the three-body force. In Sec.|nT]we discuss the 
inclusion of the Coulomb interaction using the method of 
screening and renormalization and demonstrate its valid- 
ity numerically. Section IIVI presents some characteristic 
effects of Coulomb and Urbana IX 37V force in pd elastic 
scattering and breakup. Section W\ gives our summary. 
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II. THREE-PARTICLE SCATTERING 
EQUATIONS INCLUDING IRREDUCIBLE 
THREE-BODY FORCE 

For the description of three-particle scattering interact- 
ing via three pairwise potentials v a , ex. = 1, 2, 3, we used 
Alt, Grassberger, and Sandhas (AGS) equations p~5[. In 
this section we give a short derivation of the AGS equa- 
tions including a three-body force 



(1) 



which is decomposed into three terms u a that are sym- 
metric in the exchange of particles [3 ^ a and 7 ^ a and 
can be transformed into one another by a cyclic permu- 
tation; these properties will be used later for the sym- 
mctrization of equations. The full resolvent is 

G = (E + iO-H ~J2v 7 -J2u J )- 1 , (2a) 

7 7 

G = Gp + Gp X)(W7 + %) G > (2b) 
7 

where 5p a = 1 — Sp a , E is the available three-particle 
energy and Hq the three-particle kinetic energy opera- 
tor, both in the center-of-mass (cm.) system, and the 
channel resolvents are 



G 7 = (E + i0-H -v^y 1 , 
G7 = Go + GqTjGq. 



(3a) 
(3b) 



Here Go = (E + iO — Hq) _1 is the free resolvent and 

T 7 = Vy + V-yGoT-y (4) 

the two-particle transition matrix. 

The multichannel three-particle transition operators 
U/3a are defined by the decomposition of the full resol- 
vent into channel resolvents according to 



G — 5f3aG a + GpUp a G 



(5) 



Inserting Eq. ((5]) with j3 = 7 into Eq. (|2bj) . comparing it 
back to Eq. (J3J) and taking into account Eqs. (f3"a[) and ^ 
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we obtain integral equation for the multichannel three- 
particle transition operators 

(6) 

+ ^ U7 Go(l + T 7 Go)[/ 7Q . 

7 

The on-shell matrix elements {4>p\Up a \^) a ) a re ampli- 
tudes (up to a factor) for elastic ((3 = a), rearrange- 
ment (0 / /3 / a) scattering, and breakup ((3 = 0). 
The channel states \<fi a ) are the eigenstates of the corre- 
sponding channel Hamiltonian H a — Hq + v a with the 
energy eigenvalue E while \<po) describes the free relative 
motion of three particles; the dependence on the Jacobi 
momenta and discrete quantum numbers is suppressed in 
our notation. 

For nucleon-deuteron scattering it is convenient 
to consider nucleons as identical particles in the 
isospin formalism where the symmetrized amplitudes are 
Sa( < ^/3l^/3al'/ , «)- The symmetrized transition operator 
for elastic scattering is the solution of the symmetrized 
AGS integral equation 

U = PG 1 + (1 + P)u + PTGoU 
+ {l + P)uG {l+TG )U. 

We omit the spectator index that is not needed anymore. 
The basis states are antisymmetric in the pair only, the 
full antisymmetry is ensured by P — P12P23 + P13P23 
where P a p is the permutation operator of particles a 
and p. The breakup operator is then obtained from the 
quadrature 

U = (1 + P)[G^ +u + TG Q U + uG (l + TG Q )U\. 

(7b) 

For the practical solution it is convenient to introduce 
the transition operator 

X = G„ 1 +u + TG a U + uG {l + TG Q )U (8) 

such that 

U= [P + uG Q {l + P)]X, (9a) 

U Q = (1 + P)X. (9b) 

Operator X is obtained from the integral equation 

X = G 1 + TG Q PX + (1 + TG )uG Q (l + P)X. (10) 

The practical advantage of this equation is that all terms 
have the same structure with respect to the relative mo- 
mentum of the interacting pair in the final state; this 
is convenient for the interpolation that is needed to cal- 
culate PX at each iteration step. The handling of the 
3N force is discussed in the Appendix [A] Otherwise, the 
numerical technique for solving Eq. (|10[) in momentum- 
space partial-wave representation and calculating on- 
shell elements in Eqs. @ is taken over from Refs. fla.ll7j|. 



There, two equally reliable interpolation methods for the 
two-nucleon transition matrices were used. However, in 
the presence of the irreducible 3N force we use the one 
based on splines that is more convenient than the one 
based on the Chebyshev expansion. 

Finally we note that in Refs. [H, 0] nucleon-deuteron 
scattering equations of similar form were solved. The 
relation between our transition operator X and the op- 
erator T of Refs. [l|, 0| , not to be confused with the two- 
nucleon transition matrix T in the equations above, reads 
X = Ga l +T. 



III. INCLUSION OF THE COULOMB FORCE 

In order to include the Coulomb interaction we use 
the method of screening and renormalization [ID, [3] as 
described in detail in Refs. [HI, [H, for pairwise in- 
teractions. The screened Coulomb potential w 7 r, that 
in the configuration space representation has the form 
w 1 a{r) — (a e /r) exp(— (r/R) n ), a e « 1/137 being the 
fine structure constant, is added to the hadronic proton- 
proton (pp) potential w 7 . The pp transition matrix ((4]) is 
calculated with the full interaction v~ t + w 7 r and used to 
solve the AGS equations yielding multichannel transition 
operators (U^ and in the symmetrized ver- 

sion); their dependence on the Coulomb scree ning radius 
R is indicated. Following the strategy of Refs. [ill fill Il9l] 
it is straightforward to decompose the amplitudes into 
long-range and Coulomb-distorted short-range parts also 
when the 37V force is present. The only difference is 
that the integral equations for the reduced short-range 
operators U$ of Refs. [H H 03 contain additional 
three-body force terms, i.e., 

V/fa = ^Pa^aR +V a ) +U a + 5p a W a R 

+ X^W^ + u i + 8frWpR)GiRU$, (11a) 
7 

U^T = G- R + v a +u a + 5> 7 + Ul )G lR U$, (lib) 

7 

where G a R and W a R are defined in Refs. [HI, [H, Il9j . 
Nevertheless, the relation of Ujfa to the full AGS opera- 
tors lli^} as well as the screened Coulomb contributions 

pa 

diverging in the R — > 00 limit remain the same. Thus, 
the renormalization prescription can be taken over from 
Refs. 0, d d where the amplitudes for pd elastic scat- 
tering and breakup referring to unscreened Coulomb are 
calculated as 

{4>'\uM\</>) = {<f>'\Tg°-\4>) 

+ lun {Z-\ qi )(cP>\[U (R) -T^]\<t>)}, 

K—>oo 

(12a) 

{MUP\4>) = lim {z R Hpf){MU { R) \^)Z R H qi )}. 

H — >oo 

(12b) 
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The long-range part of the screened elastic scattering 
amplitude is given by the two-body on-shcll transition 
matrix {4>'\T^ m '\4>) derived from the screened Coulomb 
potential between proton and the cm. of the deuteron. 
Renormalized by Z^ i 1 (qi) in the R — ► oo limit it con- 
verges (in general, as a distribution) to the well-known 
pure Coulomb amplitude ((/>' '|T£ m \<f>) . The Coulomb- 
distorted short-range parts (</>'|[?7 (fl) - T£ m M and 

(<Po\Uq \4>) are calculated numerically at finite R since 
after the renormalization they rapidly converge with R 
due to their short-range nature; we only have to make 
sure that R is large enough for the desired accuracy. In 
Eqs. (112|) qi is the magnitude of the initial relative pd 
momentum and pj is the magnitude of the final relative 
pp momentum. The renormalization factors Zn(qj) and 
Zn(pf) are diverging phase factors defined in Ref. [l3| for 
a general screening and given in Refs. [HI, H, for the 
form of screening used there as well as in this work. We 
get well- converged results for the observables when the 
Coulomb-distorted short-range part of the amplitudes is 
calculated using the screening function with n = 4 and 
R = 10 fm (30 fm) for pd elastic scattering (breakup). 
With those values of R the partial-wave expansion con- 
verges slower than in the nd case but the AGS equations 
still can be solved in the partial-wave basis. In our first 
calculations with Coulomb [HI, [H| high two-baryon par- 
tial waves needed for the convergence were included using 
the perturbative approach of Ref. 20] that was confirmed 
to be highly reliable in later calculations where we found 
an efficient method 21] to include the high partial waves 
exactly. This latter method [2l[ is used in the present 
work where we obtain fully converged results by taking 
into account the hadronic interaction in two-nucleon par- 
tial waves with total pair angular momentum / < 5 and 
the screened Coulomb interaction in two-proton partial 
waves with pair orbital angular momentum L < 14. The 
partial waves with total 3N angular momentum J < r§- 
are considered but it is fully sufficient to include the 3N 
force only in those with J < Of course, both total 
3iV isospin T = \ and | states are included. 

Thus, the presence of the 3iV force does not change 
the rate of the i?-convergence for the observables whose 
detailed study was given in Refs. [HI, [l8[ and will not 
be repeated here. However, a recent work [22j on an al- 
ternative Coulomb treatment in pd scattering proposed 
a renormalization prescription that is different from ours 
given in Eqs. (fT2|) . According to Ref. [HJ, the pd elastic 
scattering amplitude calculated with screened Coulomb 
does not need renormalization at all, i.e., the limit 
lim/f^oo ((f>'\U( R ' \4>) should exist. The renormalization 
for the breakup amplitude is needed but it is different 
from ours given in Eq. (|12bj) . However, numerical re- 
sults of Ref. [22| involve approximations that are not 
well under control: the screened Coulomb transition ma- 
trix is approximated by the screened Coulomb poten- 
tial and, furthermore, particular contributions to the pd 
scattering amplitudes that include first order terms in 
the screened Coulomb transition matrix are neglected. 



Therefore we feel the need to clarify the issue of renor- 
malization. We study the dependence on the screening 
radius R for selected components of the nonrenormal- 
ized and renormalized pd elastic scattering amplitudes, 
{(f)'\U {R *>\(j)) and (ct)'\U^\(j}), calculated at proton lab en- 
ergy E p = 9 MeV; the AV18 NN potential [|[ and the 
Urbana IX 3N force [§] are taken as the hadronic in- 
teraction. In Fig. [1] we show spin-nondiagonal ampli- 
tudes that have only the Coulomb-distorted short-range 
part. The nonrenormalized amplitude (<f)'\U( R '\(f>) shows 
a clear independence; in fact, its absolute value becomes 
i?-independent but the phase is proportional to In R. In 
contrast, the renormalized amplitude {4>'\U^\(p), within 
the accuracy of the plot, becomes independent of R for 
R > 10 fm. The renormalized spin-diagonal amplitude 
that has also a long-range part becomes i?-independent 
as well, as Fig. [2] demonstrates. Thus, fully converged 
numerical results without uncontrolled approximations 
clearly support the standard screening and renormaliza- 
tion theory as given in Eqs. (Ti"2"|) and not the one of 
Ref. [2|. 

IV. RESULTS 

Numerical results of this paper are derived from the 
AV18 NN potential [|[ together with the Urbana IX 3N 
force @; this combination, AV18 + UIX, is one of the 
most widely used NN + 3N force models. Hadronic 
charge dependence is fully included, i.e., pp and np po- 
tentials are used within the isospin formalism. Point 
Coulomb interaction is added for pp but the additional 
electromagnetic terms of AV18 are not included. In or- 
der to isolate Coulomb and 3N force effects, calculations 
without Coulomb or without the 3N force are performed 
as well. Obviously, we have many more predictions than 
it is possible and wise to show. Therefore we make a ju- 
dicious selection and present the most interesting cases. 
The readers, dissatisfied with our choice, are welcome to 
obtain the results for their favorite data from us. 

In Fig. [3] we show results for the differential cross 
section and analyzing powers of pd elastic scattering at 
9 MeV proton lab energy (18 MeV deuteron lab energy). 
The Coulomb effect is important in the whole kinematical 
regime. There are well-known discrepancies @ around 
the maximum of the vector analyzing powers A y (N) and 
iTn and the first minimum of T%y. The 3N force effect 
is rather small; it is mostly a scaling effect due to the 
change in the binding energies of 3 H and 3 He that are 
calculated in the Appendix [Bl In addition, observables 
shown in Fig. [3] serve as a benchmark since they already 
have been calculated in Ref. [6( using AV18 and AV18 
+ UIX force models together with the Coulomb force. 
The 3N Schrodinger equation was solved in configuration 
space using the Kohn variational principle and explicitly 
imposing the proper Coulomb boundary conditions. The 
pd elastic scattering results obtained with the two meth- 
ods for including the pp Coulomb interaction, i.e., our 
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FIG. 1. (Color online) Real and imaginary parts of nonrenormalized (left side) and renormalized (right side) pd elastic scattering 
amplitudes (in arbitrary units) at E p = 9 MeV shown as functions of the cm. scattering angle. The initial and final particle 
spin projection quantum numbers are m p = |, m' p = — |, and mj = m' d — 1. Results for the AV18 + UIX force model obtained 
with the screening radius R = 5 f m (dashed curves), 10 mi (dash-dotted curves), 15 fm (dotted curves), and 20 fm (solid curves) 
are compared. 



momentum-space screening and renormalization method 
and the one of Ref. [f| , were compared in Ref. [24[ tak- 
ing the AV18 NN potential as the hadronic interaction 
alone; a good agreement over a wide range of energies 
for all studied observables was found. Comparing Fig. [3] 
to the 9 MeV results of Ref. [6( it is easy to see that 
the predicted 3iV force effect is the same in both cases. 
Thus, the agreement between the momentum-space and 
configuration-space results including both the 37V force 
and Coulomb and is as good as without the 3iV force. 

In Fig. [4] we show the differential cross section and an- 
alyzing powers for pd elastic scattering at 135 MeV pro- 
ton lab energy (270 MeV deuteron lab energy) where the 
configuration-space calculations are not available so far. 
As found in Ref. [ll[ and confirmed in the present work, 
the Coulomb effect at this relatively high energy is large 
only at forward angles. We therefore do not show it sep- 
arately. Instead, since the 3N force effect is quite signif- 
icant, it is interesting to compare two different 3A^ force 
models: (a) Urbana IX, based on the Fujita-Miyazawa 
force {27), i.e., the two-pion (27r) exchange with an in- 
termediate A-isobar excitation, and supplemented by a 
purely phenomenological repulsive short-range part; (b) 
an effective 3N force due to explicit A-isobar excitation 
that uses no static approximation for the propagation 
of the A, includes beside the pion also the exchange of 



heavier mesons p, ui, and a, and has higher order con- 
tributions, e.g., three-meson ring diagrams; all contribu- 
tions are consistent with each other since they are built 
from the same two-baryon coupled-channel potential CD 
Bonn + A 3 that is as realistic as its purely nucleonic 
reference potential CD Bonn (28|. The A-isobar effect 
is isolated as the difference between predictions of the 
CD Bonn + A and CD Bonn potentials. In most cases 
both Urbana IX and A-isobar effects are qualitatively 
similar although the former is usually larger, especially 
around the minimum of the deuteron vector analyzing 
power Ay(d) where the predictions of AV18 + UIX are 
considerably closer to the experimental data. The reason 
why the A-isobar effect is smaller, at least partially, is the 
NN dispersion that is generated by the explicit A-isobar 
excitation and often competes with the effective 3N force 
Q. The two 3N force models show qualitatively differ- 
ent behavior around the first minimum of the deuteron 
tensor analyzing power A xx where the CD Bonn + A 
accounts for the data better. Finally we note that our 
AV18 + UIX results without Coulomb (not shown here) 
are consistent with the corresponding results of Ref. [l| ■ 
In Fig. \5\ we show the fivefold differential cross sec- 
tion for pd breakup at 13 MeV proton lab energy in 
few special kinematical configurations that are charac- 
terized in a standard way by the final-state polar angles 
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FIG. 2. (Color online) Real and imaginary parts of the renor- 
malized pd elastic scattering amplitude. The initial and final 
particle spin projection quantum numbers are m p = m' p — | , 
and m,d = m' d = 1. Curves as in Fig. [1] 



of the two detected protons and by the azimuthal angle 
between them, (0i,02,yi2 — V>2 — <fi)- Though the in- 
clusion of Coulomb slightly improves the agreement with 
data in the space star configuration, the Coulomb ef- 
fect is far too small to reproduce the difference between 
the experimental pd and nd data and to resolve the so- 
called space star anomaly. Slightly larger and beneficial 
Coulomb effects are seen in quasifree scattering (QFS) 
and collinear configurations; the differential cross section 
is decreased around QFS peak and increased around the 
collinear point (minimum) and np final-state interaction 
(FSI) peaks. The remaining discrepancies around the 
7ip-FSf peaks, at least to some extent, may be due to 
the finite geometry, not taken into account in our calcu- 
lations owing to the lack of information on experimental 
details. The 37V force effect is very small except for the 
QFS configuration. All these findings are consistent with 
our previous results [l8| derived from the CD Bonn + A 
and CD Bonn potentials. 

Finally, we consider pd breakup at 130 MeV deuteron 
lab energy that was measured recently in a variety of 
kinematical configurations [12, [H| . In some of them we 
found sizable Coulomb effects for the differential cross 
section [18j, |33| but very small A-isobar effects. In con- 
trast, the present calculations as well as those without 
Coulomb given in Ref. [32j reveal visible effects of the 
Urbana IX 37V force. Therefore a more extensive study 



of the interplay between the Coulomb and the 37V force 
is needed. In Fig. [5] we show the fivefold differential cross 
section for few kinematical configurations with small rel- 
ative azimuthal angle (fi 2 - The central point of the 
(13°, 13°, 20°) configuration corresponds to very low rel- 
ative pp energy in the final state, E pp < 0.2 MeV, where 
the differential cross section and also the 37V force ef- 
fect are strongly reduced by the Coulomb as a result of 
the pp repulsion; the found Coulomb effect is well sup- 
ported by the experimental data. The relative pp en- 
ergy in the final state increases with the relative polar 
and azimuthal angles of the protons and therefore the 
Coulomb effect decreases in the remaining configurations 
of Fig. [5] but the 37V force effect becomes more visi- 
ble. The relative pp energy gets larger at <p\2 = 160° 
in Fig. [7] leading to an increase of the differential cross 
section due to the Coulomb, especially in the configu- 
rations with smaller Since the total breakup cross 
section at this energy is almost unchanged by Coulomb 
as demonstrated in Ref. [l8| and confirmed in the present 
work, one may expect in particular configurations an in- 
crease of the differential cross section due to Coulomb to 
compensate for the strong decrease in the regions with 
low Epp. The inclusion of the Urbana IX 37V force, as in 
Fig. O increases the differential cross section, especially 
in the configurations where both 9i are large. Thus, at 
smaller 9i Coulomb is dominating while at larger 9i the 
Urbana IX 37V force becomes more significant. In some 
kinematical regimes, e.g., around the central peaks of 
the configurations 20° < 0i < 30°, 6 2 = 20°, tp 12 = 160°, 
both are equally important leading to quite a satisfactory 
description of the data in all studied configurations. 

Deuteron analyzing powers were measured in the same 
experiments [32], [331 ] . however, the data analysis is not 
yet completed. We therefore present only few examples 
in Fig. [8] demonstrating that for the spin observables like 
A xx Coulomb and 37V force effects may take place in 
completely different regions of the phase space compared 
to the differential cross section. E.g., A xx remains un- 
affected by the Coulomb in the (13°, 13°, 20°) configura- 
tion with very low E pp where a strong decrease of the 
cross section was found in Fig. [51 but shows a moderate 
Coulomb effect in the (25°, 20°, 40°) configuration where 
it was almost negligible in the case of cross section. In 
contrast, the 37V force effect in the latter configuration is 
only significant for the differential cross section but not 
for A xx . The last configuration in Fig. [8] shows moderate 
effects of both Coulomb and 37V force competing with 
each other unlike in the case of the differential cross sec- 
tion in Fig. \7\ A detailed study of deuteron analyzing 
powers in more kinematical configurations is postponed 
till the finalization of the experimental data that is ex- 
pected soon [34j]. 
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FIG. 3. (Color online) Differential cross section and analyzing powers for pd elastic scattering at 9 MeV proton lab energy as 
functions of the cm. scattering angle. Results for AV18 (dashed-dotted curves) and AV18 + UIX (solid curves) force models, 
both including Coulomb, and for AV18 + UIX without Coulomb (dotted curves) are compared with the experimental data 
from Ref. 



V. SUMMARY 

In this paper we derive AGS integral equations 
for proton-deuteron scattering including an irreducible 
three-nucleon force and solve them in the momentum- 
space partial wave representation. The Coulomb interac- 
tion between the protons is included using the screening 
and renormalization method whose validity is confirmed 
by the numerical results for scattering amplitudes; they 
are well converged with respect to the screening and with 
respect to the quantum number cutoffs. AV18 NN po- 
tential with the Urbana IX 3iV force are used as the 
hadronic interaction model to calculate the observables 
of elastic scattering and breakup. In the low-energy pd 
elastic scattering, where configuration-space calculations 
Q with the same dynamic input are available, a good 
agreement between our results and those of Ref. [a] is 
found. For pd elastic scattering at higher energies and for 
pd breakup we provide first results for the AV18 + UIX 
force model with Coulomb. In higher-energy pd elastic 
scattering where the Coulomb effect is confined to small 
scattering angles but the 3N force effect is significant, we 
compared two models, Urbana IX and the effective 3N 
force due to explicit A-isobar excitation, and with few 
exceptions found a qualitative agreement between them. 



In breakup the Coulomb effects are full y co nsistent with 
those found in our earlier calculations [ill . [l8j with dif- 
ferent hadronic interactions. The inclusion of Coulomb is 
unable to resolve the space star anomaly at low-energies, 
but clearly improves the description of the experimen- 
tal data at 130 MeV deuteron lab energy where, de- 
pending on the kinematical configuration, it may signifi- 
cantly decrease or increase the differential cross section. 
A moderate Urbana IX 3N force effect is seen as well 
that increases the cross section which is only slightly un- 
derpredicted by the theory. A complicated interplay of 
Coulomb and 3N force effects takes place in the spin ob- 
servables; they will be studied extensively in the future 
using other 3N force models as well. 
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FIG. 4. (Color online) Differential cross section and analyzing powers for pd elastic scattering at 135 MeV proton lab energy. 
Results for AV18 (dashed-dotted curves), AV18 + UIX (solid curves), CD Bonn (dotted curves), and CD Bonn + A (dashed 
curves) force models, all including Coulomb, are compared with the experimental data from Ref. [25| (full circles) and Ref. [26l ] 
(open squares). 




FIG. 5. (Color online) Differential cross section for pd breakup at 13 MeV proton lab energy in space star (left), quasifree 
scattering (middle), and collinear (right) configurations as function of the arclength S along the kinematical curve. Curves 
as in Fig. [3] and the experimental data from Ref. [29( (full circles). For the space star configuration also the nd data from 
Refs. |30l . 1311 ] (open and full squares) are shown. 



Appendix A: Urbana IX 3iV force 

The Urbana IX 3N force has 27r-exchange and phe- 
nomenological repulsive short-range terms, 

a/37 cyclic 

+ C2„[X*p, Xp 7 ][T a •Tp,Tp-T J ] + UqT^T^) , 

(Al) 



where curly and square brackets denote anticommutator 
and commutator, respectively. The strength constants 
are A 27I = -0.0293 MeV, C 27r = jA 27n and U = 0.0048 
MeV; the latter should be not confused with the sym- 
metrized breakup operator Uq defined in Sec. [TTJ r Q is 
the isospin vector operator of the nucleon a, and X*g 
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FIG. 6. (Color online) Differential cross section for pd breakup at 130 MeV deuteron lab energy in selected kinematical 
configurations with small relative azimuthal angle. Curves as in Fig. [3] and the experimental data from Refs. [H, |3^]. 




S (MeV) S (MeV) S (MeV) 

FIG. 7. (Color online) Differential cross section for pd breakup at 130 MeV deuteron lab energy in selected kinematical 
configurations with large relative azimuthal angle. Curves as in Fig. [3] and the experimental data from Ref. [3^ . 

and T a p are local potential-like two-nucleon operators de- of X™ g and T^g to the momentum space using spheri- 

pending on the relative coordinate of the nucleons a and cal Bessel functions is straightforward. For decomposing 
(3 and, in the case of X£ g , also on their spins; the explicit 
expressions can be found in Ref. [9( . The transformation 
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FIG. 8. (Color online) Deuteron analyzing power A xx for pd breakup at 130 MeV deuteron lab energy in selected kinematical 
configurations. Curves as in Fig. [3] 



(|Aip into three symmetric parts we follow Ref. [35|, i.e., 

u 7 = A2nX%p{X^[2T a ■ r 7 - % -T a ■ (t/3 x r 7 )] 

+ X; a [2r 7 • r p + l -r a ■ (t x t 7 )] } (A2) 

a/?7 being cyclic. In the 37V scattering and bound state 
equations, (fT0|) and (|B3|) . the 3N force component (|A2[) 
always acts on a state |y) that is fully antisymmetric 
due to (1 + P). Furthermore, taking into account that 
X£g, T a /3, and u 7 are symmetric under exchange of the 
nucleons a and /?, one can prove that certain terms in 
Eq. (|A2j) yield equal contributions [35[ and calculate the 
desired matrix elements of the 3/V force component as 

^1 ^2 ^3 

x 7 {fi\[2r a ■ r 7 - -r a ■ (r fj x t 7 )]|j/ 2 )q 

1/1 1/ 2 1^3 

X a (v2\Tp^\l>3)aa(va\y), 

(A3) 

where each of the intermediate states \vi) a , with Vi ab- 
breviating all continuous and discrete quantum numbers 
and a denoting the spectator, is antisymmetric under 
the exchange of the nucleons and 7. The transfor- 
mation from the basis with the spectator a to the ba- 
sis with the spectator 7 is needed in order to evalu- 
ate all the matrix elements of X*p and in their 
proper bases. Since \y) is fully antisymmetric, a (u3\y) = 
i(v3\y)- Furthermore, 1 {vi\v 2 ) a = \-,{vi \P\v 2 ) 1: and 
7 {vi\[2T a ■ r 7 - ^r a ■ (173 x T 1 )]\v 2 ) a differs from 7 (^i|^) a 
only by the isospin factor that can be found in Ref. [35j |. 
Thus, the calculation of the 27r-exchange and short-range 



terms needs only one basis transformation for each of 
them, and all intermediate states are antisymmetric with 
respect to the pair. This constitutes the advantage of 
the above method [35[ for handling the Urbana IX 3A?" 
force over the one used in Refs. [l| that needed two ba- 
sis transformations for each term involving, in addition, 
intermediate nonphysical symmetric states. Finally we 
note that for handling other types of 3N forces to be 
used in the future calculations the technique proposed in 
Ref. [36[ seems to be very promising. 

Appendix B: Three- nucleon bound state 

Starting from the Schrodinger equation in the integral 
form, 

l*> =GoX>« + OI*>> (Bl) 

a 

decomposing the 3N bound-state wave function \^>) = 
S« IV'q) m t° hs Faddeev components \ip a ) with a — 
1, 2, 3, and using the definition of the two-particle transi- 
tion matrix (J4| it is straightforward to derive the Faddeev 
equation 

\rl> a ) = G T a 2 S Pa \i> p ) + (1 + G T a )G u a ^ |^). 



(B2) 

The symmetrized version of it reads 

|V> = G TP\i,) + (1 + G T)G u(l + P)\i>) (B3) 

where the bound-state wave function is obtained as \^>) = 
(1 + P)\ip}- In practical calculations we solve the Fad- 
deev equation for \x) = Gq 1 ]^) that has an advantage of 
having exactly the same form of kernel as the scattering 
equation (TTU]). In contrast to scattering calculations, we 
include hadronic interaction in two-nucleon partial waves 
with total pair angular momentum / < 6 and take into 
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account all electromagnetic terms of the AV18 NN po- 
tential. Our results for binding energies and wave func- 
tion probabilities of 3 H and 3 He nuclei with and without 
Urbana IX 3iV force are collected in Table Q] they are in 
good agreement with the results of other groups [37], [38| . 



\E B \ (H ) P s , P P P c 



T=3/2 



AV18 7.621 46.72 1.292 0.066 8.509 0.0025 

AV18 [37] (r) 7.624 46.73 1.293 0.066 8.510 0.0025 

AV18 [37] (p) 7.621 46.73 1.291 0.066 8.510 0.0025 

AV18+UIX 8.478 51.28 1.055 0.135 9.302 0.0025 

AV18+UIX [37](r) 8.479 51.28 1.054 0.135 9.301 0.0025 

AV18+UIX [37j(p) 8.476 51.28 1.052 0.135 9.302 0.0025 

AV18 6.923 45.68 1.526 0.065 8.466 0.0081 

AV18 [37] (r) 6.925 45.69 1.530 0.065 8.467 0.0081 

AV18 [37] (p) 6.923 45.68 1.524 0.065 8.466 0.0081 

AV18+UIX 7.748 50.21 1.239 0.132 9.249 0.0075 

AV18+UIX [3£] (r) 7.750 50.21 1.242 0.132 9.248 0.0075 

AV18+UIX [3£| (p) 7.746 50.21 1.235 0.132 9.248 0.0075 



TABLE I. Absolute value of binding energy, expectation 
value of kinetic energy (both in MeV), and probabilities of 
the wave function components (all in %) for 3 H (top) and 
3 He (bottom) nuclei calculated with AV18 and AV18+UIX 
force models. Results of Ref. (37| obtained using coordinate- 
space (r) and momentum-space (p) frameworks are listed as 
well. 
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